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Abstract 

We give an implicit equation for the accessory parameter on the torus which is 
the necessary and sufficient condition to obtain the monodromy of the conformal 
factor. It is shown that the perturbative series for the accessory parameter in 
the coupling constant converges in a finite disk and give a rigorous lower bound 
for the radius of convergence. We work out explicitly the perturbative result to 
second order in the coupling for the accessory parameter and to third order for 
the one-point function. Modular invariance is discussed and exploited. At the non 
perturbative level it is shown that the accessory parameter is a continuous function 
of the coupling in the whole physical region and that it is analytic except at most a 
finite number of points. We also prove that the accessory parameter as a function 
of the modulus of the torus is continuous and real-analytic except at most for a zero 
measure set. Three soluble cases in which the solution can be expressed in terms of 
hypergeometric functions are explicitly treated. 



1 Introduction 



Liouville theory plays an important role in several fields both at the classical and quantum 
level [H|2j[3j|ll|5l[6l[7J[8j[9]. R- ecen ^y a renewed interest has developed due to a conjecture 
[TOT ITT] that Liouville theory on a Riemann surface of genus g is related to a certain class of 
N = 2, 4-dimensional gauge theories and the conjecture has been supported by extensive 
tests on genera and 1 pJJ [121 H3] and proven in a class of cases [HI [15] . At the classical 
level the key point in solving the theory is the determination of the accessory parameters 
which on the sphere are related to the semiclassical limit of the operator product expansion 
via the Polyakov relation. 

The determination of the accessory parameters turns out to be a highly transcenden- 
tal problem. The mathematical literature is concentrated mainly on the limit case of 
parabolic singularities i.e. punctures. On the other hand in quantum Liouville theory, 
elliptic singularities which exhibit a continuum spectrum are of most interest. 
In the three point problem on the sphere the accessory parameters are algebraically fixed 
by the Fuchs relations. On the other hand the four point problem on the sphere [H 
[TBI [TTl [T8] and the one point problem on the torus [191 EDI EH [15] lead to differential 
equations with four regular singularities which are special cases of the Heun equation. 
Higher number of point lead to still more complex equations. 

The Heun accessory parameter depends on three quantities: the coupling rj whose 
physical range is < i] < 1/2, the modulus r and a scale parameter. The dependence 
of (3 on the scale parameter is trivial while (3 turns out to be a weight two modular form 
with some simple conjugation and inversion properties. This allows to predict through 
invariance argument the value of /3 in two special cases [191 EH1 EJJ: 1) The so called 
harmonic case i.e. the square; 2) The equianharmonic case i.e. the rhombus with opening 
angle tt/6. In both cases the value of is zero. In these two cases the Heun equation 
reduces through respectively a quadratic and a cubic transformation to an hypergeometric 
equation and thus the conformal factor can be explicitly given in terms of hypergeometric 
functions for any value of the source strength in the physical region [201 EI]. Such a 
reduction is possible due to a symmetry in the parameters ei,e2,e3 which together with 
infinity give the position of the singularities in a two sheet plane which describes the torus. 
In order to connect the Heun equation to more familiar cases Maier [22] examined all 
rational substitution of the independent variable which transform the Heun equation into 
an hypergeometric equation. He found that this transformation occurs only for certain 
polynomial of degree 2,3,4,5,6 . The harmonic case corresponds to the order two and 
the equianharmonic case to an order three case. The other transformations introduce in 
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addition to the physical source some additional "kinematical sources" which correspond 
to spurious sources. Thus despite the interest of the transformation, as far as the single 
source problem on the torus is concerned, no new physically interesting case is reached. 
The nature of the dependence of the accessory parameters on the moduli and the source 
strengths is not completely known. The reason is that while the proof of the uniqueness of 
the solution is relatively simple [23J 121] the existence of the solution relies on a variational 
method i.e. on the minimization of a certain functional |2"E]. In [22] this was achieved 
by expanding the conformal factor in terms of a complete set of functions while the 
more modern treatment of [26J exploits the techniques of Sobolev spaces proving first the 
existence of a weak solution and then the existence of the solution. The outcome is that 
it is very difficult at the end to follow the nature of the dependence of the solution on the 
coupling and the moduli. 

An exception is the case of parabolic singularities (punctures) were general properties of 
fuchsian mappings can be applied. Using such a technique Keen, Rauch and Vaquez [T9"] 
found that the accessory parameter for the torus with one parabolic singularity (puncture) 
is a real-analytic functions of the modulus; in addition in [19] some numerical investigation 
of the accessory parameter was performed. Zograf and Takhtajan [27] treated the case 
of parabolic singularities on a Riemann surface of genus 0. The result of [191 12Z] is that 
the accessory parameters are real-analytic functions of the moduli. Kra [28J gave an 
extension of such a result to the case of a collection of parabolic and a special class of 
elliptic singularities i.e. finite order elliptic singularities where the strength of the source 
can assume only the values rj = 1/2(1 — 1/n), n G Z + . This is a discrete set which 
accumulates to the parabolic point. On the other hand in quantum Liouville theory, 
elliptic singularities which exhibit a continuum spectrum are of most interest. 
In the general case of elliptic singularities and parabolic singularities it was proved in 
[29~1 [30] that the accessory parameters are real-analytic in the couplings and in the moduli 
in an everywhere dense open set: given a value of the coupling, if the accessory parameter 
is not analytic at that point there is an open set as near as we like to the given point, on 
which the accessory parameter is analytic. 

Here we shall prove a much stronger result i.e. that the accessory parameter for the 
torus is an analytic function of the coupling in the whole physical region except at most 
a finite number of points and it is a real-analytic function of the modulus in the whole 
fundamental region except a zero measure set. 

In proving such results we shall rely on some properties of the solution which are extracted 
using potential-theory techniques which were used in the solution of the uniformization 
problem, combined with some results on analytic varieties [3T]. The first is the existence 
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and uniqueness property of the solution, a result which goes back to Picard himself. The 
second is the boundedness property of the solution (f> of the Liouville equation and its first 
and second derivatives with respect to the argument, in any region which excludes finite 
disks around the singularities and which was proved in [30] . 

The accessory parameter /3 obeys an implicit equation. From such implicit equation 
a power series expansion for /3 in the coupling 77 can be extracted and we prove such 
expansion to be rigorously convergent in a finite disk. We also compute a rigorous lower 
bound on the convergence radius. We compute also explicitly the expansion of in 77 up 
to second order in terms of integrals of elliptic and related functions. 
For general couplings i.e. couplings not necessarily small, exploiting the uniqueness the- 
orem and some results on complex analytic varieties [31J we are able to prove that the 
which solves the monodromy problem is analytic in the whole physical range of the 
coupling except at most for a finite number of points. 

The nature of the dependence of the accessory parameters on the moduli of a punctured 
Riemann surface is important in several respects; e.g. the C l nature of such a dependence 
is an essential input in proving Polyakov relation on the sphere [291 USB EH E2] • Here we 
prove that both for elliptic and parabolic singularity the dependence of (3 on the modulus 
is real-analytic except for a zero measure set thus extending the results of [291 E0] ■ The 
technique developed here can be applied to the four or higher point functions on the 
sphere and also to higher genus surfaces. 

Within the AGT [TOl ITT] correspondence Ferrari and Piatek [18] exploited the relation 
between the semiclassical limit of quantum Liouville theory and the Nekrasov-Shatashvili 
limit of the N = 2, U(2) super Yang-Mills theory to give an expression of the accessory 
parameter for the 4-point function on the sphere in terms of a contour integral containing 
the ratio of the column length of critical Young diagrams. 

It should be possible to extend such technology to the case of the torus. On the other 
hand once this is accomplished, a direct comparison with the result obtained here will 
not be straightforward as they are based on different expansions. In the present paper, 
the accessory parameter has been considered as a function of the source strength and an 
expansion in the source strength given. Instead in the approach of [18] an expansion of 
the accessory parameter in the position x of the fourth singularity w.r.t. the position of 
the first z = singularity appears. A similar approach which computes the accessory 
parameter expanding in x is found in [TJ. As for the torus the modulus is related to the 
positions of the singularities in the w-plane, it appears that the such expansion should 
correspond to a perturbation around the degenerate case in which two singularities co- 
incide i.e. the infinite strip which we treat in section [HI In [20J a general perturbation 
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technique under the variation of the moduli has been developed, and one could apply it 
to the present situation allowing a direct comparison. 

The paper is organized as follows: In section [2] we write down the differential equation 
in the cut-it plane and derive the monodromicity condition. The fulfillment of a single 
complex implicit equation is necessary and sufficient to assure the monodromic behavior at 
all singularities. In section [3] we give the explicit expression of the monodromy matrices. In 
section|l]we discuss modular invariance and the consequent determination of the accessory 
parameter in two soluble cases. Section [5] is devoted to the the relation of the action in 
different coordinate systems. Section [6] gives the exact expression of the conformal factor 
and of the one point function for three soluble cases, one of which is the limit case of 
the infinite strip. In section [7] we develop perturbation theory in the coupling constant 
up to the second order and give lower bounds for its convergence radius. In section [8] we 
give a different approach to perturbation theory by working directly with the conformal 
factor. Here we are able to go easily to third order even if the control of the convergence 
property of the series relies on the results of section [7J In section we derive the general 
analytic properties of the accessory parameter both in the coupling and in the modulus. 
In section [10] we give some concluding remark and point to some open problems. 



2 The differential equation and the monodromy con- 
ditions 

The equation 

-d z d- z( j) + e* = 2iTr ] 5 2 (z-z t ) (1) 

does not contain information about the torus. They have to be put in through periodic 
boundary conditions. 

To have a faithful representation of the torus we have to use the two-sheet representation 
of the torus in the variable u = p(z). For simplicity and without loosing generality due 
to translational invariance we shall set in this section z t = 0. 

We recall [2DI EI] that the problem of finding a solution to eq.([T]) can be reduced to 
finding the value of the accessory parameter and and of a real parameter k, such that 
the expression 

e^ /2 = [K- 2 ^tu) yi (u) - K 2 ^(u)y 2 (u)] (2) 

v2|wi 2 | 

is monodromic. yi, y 2 are two solutions of an ordinary differential equation which contains 
the parameter (3 and w\ 2 = yiy' 2 — y'\Vi is the constant Wronskian. 
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Such equation in u is given by [20] 



y" + Qy = y" + (Q + g)y = o (3) 

where 

3 ( 1 1 1 2e i 

QM = Tg ( 77. 7^ + 77. 7^ + 77. 7^ + 



16\.(^-ei) 2 (w-e 2 ) 2 (u-e 3 ) 2 (ei - e 2 )(e 3 - ei)(« - e x ) 
2^ 2e 3 



(e 2 - e 3 )(ei - e 2 )(u - e 2 ) (e 3 - e x )(e 2 - e 3 )(tt - e 3 ) 
and 

1-A 2 m + _l-A 2 it + £_tt + 

g(M) " 16 ( M - ei )( u -e 2 )( U -e 3 ) ~ 4 b'(^)] 2 = U 
where e k = p(u k ). 

We recall that A = 1 — 2rj and that rj has to satisfy the Picard condition < r\ < ~, the 
lower limit being due to the negative nature of the curvature in the bulk and the upper 
to the local finiteness of the area. The upper limit correspond to a puncture, or parabolic 
singularity. The range of e is < e < j. 
We know two independent solution to eq. 

y" + Qoy = o . (5) 

They are 

yi 0) = [4(« - ei )(« - ei )(« - ei )]V4 = = n(«), (6) 

2/f = (z - u; 3 )[4( U - ei )(« - e x )( U - ei )] 1/4 = Z b'(z)] 1 / 2 = Z U(u) . (7) 

Defining 

K(u,u') = yf\u)9{u } u') q ^-y$\u') - (u)9(u, u>) q -^±yf\u') (8) 
we can solve the full equation by the convergent expansion 

yj = (1 + K + KK + KiOC + . . . )^ (0) . (9) 

It will be useful to write 



and now the are given by 



AO) = 1 + e / (z'-z)(P + p(z'))dz'+ (11) 
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+ e 2 



(z' -z)(/3 + p(z'))dz' f (z" - z')(/3 + p{z"))dz" + ... 

LU 3 J(J 3 



and 



f 2 (z) = Z + e ( z '-z)((3 + p(z'))Z'dz'+ (12) 

Jlu 3 

+ e 2 Hz' -z){fi + p(z'))dz' f (z" - z')(/3 + p(z"))Z"dz" + ... 



Ul 3 JU) 3 



with u 3 = u>i + We shall be interested in the values of fj and their derivatives gj = /• 
at the points cj±, u 2 and cj 2 — cu±, u\ — u 2 . Due to the triangular structure of the multiple 
integrals and the boundedness of p(z) along the integration paths U3 — vuj\ and U3 — ruji we 
see that such fj and gj are analytic function of e and (3 as the series converges absolutely 
for any given e and /3. This is a well known fact. 
We have 

pUJk 

g x (u k ) = -E {P + p{z'))h{z')dz' (13) 



and 



g 2 (u k ) = 1 - e / (/3 + p{z'))f 2 {z')dz' (14) 

and similarly for the other values of the argument. The fj(oJk) and the gj{oJk) are holo- 
morphic functions of e and j3. 

There are two ways to impose monodromy. The first |20j is to exploit the symmetry for 
z t — of the equation ([T]) and of its solution under the inversion z — > — z. As proven in 
[20] it is the sufficient to impose monodromy around the three kinematical singularities 
Cj = p{ujj) under a full turn in u which corresponds to half turn in z. A more general 
method is to impose monodromy under two independent cycles of the torus [21]. While 
the first method requires the solution of the differential equation on the tracts (a; 3 ,a;i), 
(co> 3 ,Co> 2 ) the second method requires the solution on the longer tracts (co>3,u;3 — 2a;i), 
(us,U3 — 2u 2 ), i.e. along the full cycles. 

In order to derive non perturbative results we shall need some results from the potential 
theory approach to the Liouville equation. The input we shall use is Picard result [23] 
about the existence and uniqueness of the solution of the uniformization problem of which 
eq.Q is a particular case. Picard was concerned only with elliptic singularities i.e. 77 < \. 
Later the treatment was extended to elliptic and/or parabolic singularities in [2"4"| [25j . For 
a more modern treatment using Sobolev spaces see [25] . 

In [20] after choosing the canonical pair of solutions y\ it was proven in the first approach 
that the monodromies at ei, e 2 have the form 

. iidndn + bjCj) —2iajbj , 
MM=[ { ' ' / 3) 3 3 ) (15) 

ZiCjdj —t{ajdj + bjCj) 
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with a,j,bj,Cj,dj elements of a SL(2,C) matrix. We have still at our disposal a scale 
transformation on the canonical solutions under which the matrices go over to 

tiifljdj + bjCj) —2ia,jbjK~ 2 
2iCjdjK 2 —i(djdj + bjCj) 

Picard existence theorem tells us that for each j, cijbj and Cjdj are either both zero or both 
different from zero. From this remark it follows that necessary and sufficient condition 
for the existence of which renders 

djbjK^ 2 = CjdjR 2 (17) 

is that in (IT5|) 

M 12 (wi)M 21 (w 2 ) = M 12 (w 2 )M 21 (w 1 ) . (18) 

It is not difficult to prove [20J that once relation ( ITT)) is satisfied it follows that all the 
monodromies become SU(1, 1) i.e. the conformal factor tp becomes single valued and 
regular. Thus (3 has to be chosen as to satisfy (fTBl . 

Similar considerations hold in the cycle approach [21] where the necessary and sufficient 
condition for the realization of the monodromic solution takes the form 

M 12 {C X )M 21 (C 2 ) = M 12 (C 2 )M 21 (C 1 ) . (19) 




3 Computation of the monodromies 

Given the complex 

Y(u) = I 



Y(u) = I ^ | (20) 



the monodromy matrices are defined by 

Y(u) = MY(u) (21) 

where Y(u) denotes the complex after a complete turn in u at e\ or e 2 , in the first 
approach, or after a cycle, in the second approach. We have also 

Y'(u) = MY'(u) (22) 

and thus 

£i v'A (yi y'i 



Vi y 2 \V2 y 2 



M ^ (23) 
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from which 



y/2 y 2 J \-V2 yi J 

due to W12 = 2/12/2 — 2/i2/2 = 1- Application of eg. (1211) to eg. (ITT]) gives 




(24) 



M 12 (cj k ) 



2/12/1 + y[yi 



2e t2 fi(uk)9i{uk) 



(25) 



and 



M 2 i(wfc) = 2/22/2 - 2/ 2 2/2 = 2e l 2f 2 (u) k )g 2 (uJk)- 



Thus equation (ITS]) becomes in the first approach 



f 1 (u} 1 )g 1 (u 1 )f 2 (tJ2)g2(u2) = /l(W2)ffl(wj)/ 2 (wi)fli!(ll)i). 



(26) 



Using the cycle method the monodromy equation (IT9]) becomes 



#1(0)3 - 2wi)/ 2 (a;3 - 2w 2 ) = <7i(w 3 - 2u 2 )f 2 (uj 3 - 2o>i). 



(27) 



Due to the uniqueness of the Picard solution the two equations are equivalent. We shall 
come back to this property in section [71 

4 Modular invariance 

We recall some simple properties of the accessory parameter (3 which are derived from the 
differential equation and the uniqueness theorem. 
The equation 



has invariance properties related to the transformation properties of p under dilatations, 
conjugation and modular transformations [19]. From 



with n G C we have fi 2 (3(fiUi, fiu 2 ) = 0(coi,co 2 ). From p(z\ui,u 2 ) = p(z\ui,u 2 ) one 
obtains f3(ui,u 2 ) = f3(coi,co 2 ). Moreover as the lattice is left invariant under uj\ —> —uji 
and u 2 — u 2 and under oj\ -H- u 2 , p(z) is unchanged and also (3 is unchanged. Similarly 
(3 is unchanged under u\ — > coi, co 2 — > co 2 + U\. The two transformations cji — > cji, uj 2 — > 
u 2 + U\ and u\ — > u 2 , u 2 — > —oji are, apart for a dilatation, the generators T and S of the 



f( z )+e(p(z)+/3)f(z) = 



(28) 



(29) 
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modular group [12] • Thus defining (3[u2\ = (3(1, U2) we can synthesize the transformation 
properties of (3 as 

= ( C T + d) 2 (3\r} (30) 
telling us that (3 is a modular form of weight 2 and 

m=p\-T]- (si) 

From such transformation properties a few simple facts follow [T9~l |2"U| |2"T] . 

1) For Re r = 0, (3 is real; this describes the rectangle. 

2) For Re r = ±~ we have (3 = real. 

From the fact that the stabilizer in the fundamental region is ±7 except S for r = i, ST 
for r = e 2 ™ /3 and TS for r = e™ /3 [42] we have 

3) For t = i we have /3 = 0; this describes the square. 

4) For t = e 2m / 3 we also have ,5 = 0; this describes the so called equianharmonic case 
where the fundamental region is a rhombus with opening angle 2ir/6. 

5) For t = e 7 ™/ 3 we also have (3 = 0; this again describes the equianharmonic case where 
the rhombus with opening angle 27r/6 has a different orientation so it does not differ from 
case 4. 

In [20l [21] the explicit form of the conformal factor in terms of hypergeometric functions 
was given for the cases 3 and 4,5. 

From the viewpoint of the differential equation in u, modular transformations boil down 
to a simple permutation of the and a scale transformation. Thus in studying the 
monodromies in the u cut-plane with the first method we have a simple interchange of 
efc in the basic equations. If instead we exploit the cycle approach, modular invariance 
is due to the group composition properties for the transfer matrices, when we add to a 
given cycle one or more cycles. This will be relevant in discussing the modular invariance 
of the perturbation theory results. 



ar 



ct + d 



5 The action in different coordinates 

It is well known [H SI] that also the classical action has to be regularized due to the 
logarithmic divergences which arises from the kinetic term at the singularities. 
In this section we shall write the relation among the two regularized action S z and S u 
related to the z- and w-representation of the torus. 
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In the ^-representation the action is given by 

^ = J- / {\d4 A # + e*dz Adz^-^-l ft-*?- - JLr) - r? logs 2 (32) 
2ir 2tt J De 2 2 Am J dDs z - z t z - z t 

where D £ is exterior of a circle of radius e around the source at z t . Writing 

= -2<q\og\z - z t \ 2 + X + o(z - z t ) (33) 

we have also 

^± = }-f (ld(f)A d(f) + e 4 'dzAdz) 1 - -r]X + r] 2 \oge 2 (34) 
2n 2tt J Dg 2 2 

and the important relation 

^ = (35) 
2tt dr] y ' 

In order to compute the action explicitly in the soluble cases it is however better to put 

the source at the origin z t = 0, which gives rise to a singularity in u at infinity. The 

transition to the ^-representation is given by u = p{z) 

e^dz Adz = e^du A du (36) 



i.e. 



with 



Thus 



We have 



= up - log J J (37) 

J = tu = ^z)>- (38) 
= (f + log \A(u - ei) (it - e 2 ) (it - e 3 ) | (39) 



0= -2r]log|2| 2 + X + o(z) (40) 
and taking into account that w = = + o(z) we have at infinity in u 

<p=( V -^)\oguu + X^ (41) 

with 

X^ = X-log4. (42) 
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The regularized S u action now takes the form 

— S u = — / (-dip A dip + e v du A du) — 
2ix 2n J p 2 2 



1 , 3. / .du du. 1 f . du du , 

2 J Ru u u 16m J £k u-e k u- e k 



1 / 3 N r, , n r) 1 

2 ( ^2' 

1 f 1 i 

= — / (-dip A dip + e v du A du) — 
2n J D 2 ' 2 

- ^(^-^) 2 log^-(r/-^)X« +llog4-^ (43) 

where the integration in u is extended to the two sheets which describe the torus. D 
excludes disks of radius Ek around and on both sheets the exterior of a circle of radius 
R u . 

We also have 

" (44) 



27T 9?7 

Using eqs. f)34ll36|r43|) we find for the relation between the two actions 



S z = S u -\ log |( ei - e 2 )(e 2 - e 3 )(e 3 - e l )\ 2 - 2r;log2 . (45) 



consistent with eq. 



6 Soluble cases 

In this section we shall give the explicit value of the action, i.e. of the semiclassical 1-point 
function, for three soluble cases i.e. the square, the equianharmonic case i.e. a rhombus 
with opening angle 27r/6 and the limit case of the infinite strip for any coupling in the 
physical region. 

1. The harmonic case: the square 



In [201 [21] the Liouville field for the harmonic case i.e. the square was computed in terms 
of hypergeometric functions. The result was with A = 1 — 2i] 



<l>(z) 



-log^l 



K\ 



(46) 



+ loe 



1 — A1 + A3 2/w 



|k| 4 |m(z)| 



_ , 3 — A 3-)-A 5 o. . » 



n 



with 



where as usual 



From (1461) we obtain 



7 2 (i)7(^)7(H A ) 



7(x) 



r(x) 



r(i -x) 



1 — 277 

X = (3 + 477) log2-21og 7 ( — — '-) + log 7(1 -77) 
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and using (j35j) we have with 



(47) 
(48) 

(49) 



F(x) = / \og / y(x') dx' 



(50) 



1 2r? 1 

^(square) = -(3t7 + 2rf) log 2 - 4F( ' ) + F(l - 77) + 4F(~) - F(l) . (51) 



2. The 

We have [21] E 



equianharmonic case 



log(2V2|K| 



(52) 



+ lot 



1 - A 1 + A 2 3 . 
12 ' 12 ' 3' v 



\k\ a \u{z) 2[ 



,5-A 5 + A 4 



12 ' 12 '3' 



with 



IkI =9 



vrr(^)r(^) 
r 2 (|)r(^)r(i±A) 



« 3 (z)) 



(53) 



which gives 



7 4 2 77 1 — 277 77 

X = 2 log 3 + - log 2 + -77 log 2 - log 7b + {) - 21og 7 (— — ^ + log 7 (l - A (54) 
3 3 3 3 o 3 



Integrating 
S z (equianharmonic) 



2 V log 3 + h log 2 + ^ log 2 - 3F(^) + 6F(LJ3) 



- 3F(l-|)+3F(H)-6F(i)+3F(l) 



(55) 



3. The infinite strip 



1 We correct for a factor 2 in the argument of the first logarithm 
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We discuss here a limit case of the torus topology which is soluble i.e. the infinite strip. 
The vertical infinite strip is reached with the parameters e\ = 2a, e 2 = = —a. For 
a = 1 we have 



Q(u) 



1-A 2 



u + (3 



+ 



3 u 2 - 2u + 9 



(56) 



16 (u+l) 2 (tt-2) 16(m-2) 2 (m + 1) 2 ' 
The accessory parameter /3 has to be fixed to 1 otherwise the pole of order 2 at u = — 1 
would not have the correct kinematical value 1/4 as it is required for the limit of an infinite 
rectangle. It is of interest that the value 1 is already given by first order perturbation 
theory [20] 

C(w 2 )u>i - C0*>i)w 2 







(57) 



Qj 2 uJi — 001UJ2 

In fact for the case at hand i.e. e\ = 2a, e 2 = = —a the Weierstrass p and ( functions 
become 133 



u = p{z) 



-a + 



3a 



with 



and thus 



sin 2 (v / 3a z 



. r— cos(V3a z) 

(^(z) = az + v 3a — 



sin(v3a z) 
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2V3a 



lim (3 — a — 1 . 

UJ'2— >ioo 



(58) 

(59) 
(60) 



Two independent solutions of the differential equation with Q given by (|56p . canonical 
at it = 2 are 

Vl = (u - 2)y\u + I) 1 / 2 2 F l( l^, i±A, I; ^) (61) 



giving for the 

i = 
2 

with 



y 2 = (n-2) 3 / 4 (n + l) 1 / 2 2 F 1 ( 



log(3\/2) - log 



4 ' 4 '2' 3 ' 
3-A 3 + A 3 2-u, 



4 ' 2' 3 



(62) 
(63) 



1-A 1 + A 1 2 — it . 



4 ' 2' 3 
, 1 



re |u - 21 



,3- A 3 + A 3 2-u, 



2 ^V )7(i T^ : 



We find 

from which 

(strip) = - 



4 1 
X = 2r] log - + log 6 + 2 log 7(1 - 77) - 2 log7(- - 77) 



r? 2 log - + 77 log 6 - 2F(1 - 77) + 2F(- - 77) + 2F(1) - 2F(~) 



(64) 
(65) 
(66) 
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7 Second order perturbation theory and convergence 
radius 



In this section we shall develop the perturbation theory around e = 0. We shall show that 
the perturbative series is convergent in a neighborhood of e = and give a rigorous lower 
bound on the convergence radius. We shall also give the explicit expression of the first 
and second order term for the accessory parameter 0. We stress that the treatment of 
this section requires neither Picard's existence and uniqueness theorem nor other results 
from the potential theory approach to the problem. 

The perturbative series for (3 is obtained by solving the implicit equation (I26"j) or (I27|) . 
We shall use (T27D 



eF(/3, 0, e) = 9l (u 3 - 2w 1 )/a(w 3 - 2cu 2 ) - gi (u 3 - 2uj 2 )f 2 (u 3 - 2u x ) = . (67) 
We notice that 

gi (u 3 - 2u k ) = 2e{u k [3 - CM) + 0{e 2 ) (68) 
being ((z) the Weierstrass zeta- function, while 

f 2 (u 3 - 2u k ) = -2u k + 0(e) . (69) 

Thus after dividing f )27|) by e we have for the Jacobian, at e = 

d(F,F) 



J 



lQ\tU2i0i — ij02U\\~ 7^ . (70) 



e=0 



Due to the analyticity of F in /3, /3 and e we have that /3 will be an analytic function of 
e in an open neighborhood of e = [32] ■ We can therefore develop a perturbative series 
around e = 0. From (I68II69P we have [20] 

UJ 2 UJl — U)\U 2 

On equation (|71|) we can check already the following properties: 1. For uj 2 = uo\ i.e. 

the square (5\ = 0. 2. For co 2 = e 2m ^Ui i.e. the equianharmonic case /?i = 0. For 
general oji,uj 2 we verify modular invariance, i.e. invariance under the two generating 
transformations 

Ux ->■ ui, u 2 ->■ u 2 + oji (T) 
ui ->■ — u 2 , lu 2 ->■ ui (S) . (72) 

All these properties are proven easily using 

C( Wl + W2 ) = (( Wl ) + C(w 2 ). (73) 



14 



Moreover in the limit of the infinite strip we have — > 1 (see eq. (l60]) ). 
We come now to the second order. Developing eq.f lST)) we obtain for the accessory param- 
eter ft to the second order 



= A 



U>lU2 — U>\U>2 



- 2wiC(wi)C2(wi,w 3 ) - 2( 3 (u 1 ,u 3 )((u 1 ) +uj(u 1 ,u 3 ) + 



+ /3i[ - oj%uji(;(ui) + 2uxC, 3 {ujx,uj 3 ) +2utC 3 (uj l7 u 3 ) +oj 2 uf (2C(wi) +C(^2)) 

(74) 

where in the above expression A is the first order result ( I7ip and 



2 _ 2 2 - 

+ gA^fC(Wl) + g/^lWl^ + -PlP&lUi ~ {Wl <B" W 2 } 



C 2 (^ 3 ) = /*C(/)dy = log^, C 3 (z,u 3 )= I ( 2 (z',u 3 )dz' 



W3 

2 f i\ 



I(z,co 3 ) = / CV)^- (75) 
j 0)3 

Some comments are in order about such a result. £(z) is an entire function and I(z, 0*3) 
is a single-valued function of z in the fundamental parallelogram due to the absence of 
the constant term in the expansion of £(2) 

C(z) = i-^- + ... (76) 
s v ; z 2 2 ■ 3 ■ 5 V ; 

On the other hand (z(z, u 3 ) and ( 3 (z,u 3 ) are not single-valued functions; nonetheless the 
combinations 



2(3(2,^3) +zCs(z,u 3 ) (77) 
are single valued in the fundamental parallelogram. In fact we have 

£2(2, ou 3 ) = log z + regular terms (78) 

and 

( 3 (z, co 3 ) = z log z + regular terms (79) 



from which it follows that the terms £1X2(^1, W3)+C3( w i) ^3) an d ^1(3(^1^ w 3)+^iC3( w i; ^3) 
in (l74p are well defined. 

It is very important that the paths chosen in evaluating (2(^1,^3) and £3(^15 ^3) are the 
same even if there is no preferred path. 
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One can easily verify that, as expected, for the square and the equianharmonic case such 
second order contribution vanishes. On the other hand is very cumbersome to verify 
directly the modular invariance of eq. (j74p . Modular invariance of eg. f)74p is assured at the 
exact and also perturbative level by the group composition properties of the monodromies 
over cycles. Obviously the same result (1741 is obtained using the first approach to the 
monodromy problem i.e. starting from eq. (12B|) . Iterating the process in eq. (I67p one can 
go to higher orders. 

We come now to the convergence radius of the perturbative series in e. A rigorous lower 
bound on the convergence radius can be obtained applying Rouche theorem [4"5] . 
Equation flBT|) can be rewritten as 

p-p 1 + eG(pJ,e) = (80) 

with fix given by (1711) . It will be useful to exploit the polarization technique [36] i.e. to 
introduce in addition to /3 an other independent complex variable /3 C and consider the 
system 

P-Px + eGtf^ce) = (81) 
P c -Px + eG(P c ,P,e) = (82) 

where G is the analytic function obtained by conjugating in the power expansion the 
coefficients of G. Obviously if P, P c for real e is a solution of the above system, also p c , (3 
is a solution. If, always for real e, the solution is unique then we have P c = P and such 
solution is the solution of the monodromy problem. 
Given a positive constant B we can always find a 5 such that for |e| < 8 

B>\eG(P c ,P,e)\ (83) 

for all P, P c with \P — Px\ < B, \fi c — px\ < B . Then due to the analyticity of G we can 
apply Rouche theorem to conclude that (182]) for |e| < 5 has one and only one solution 
Pc = Pc(P, e), with \p c — Pi\ < B. Moreover such p c will be an analytic function of P and e. 
We substitute now such P c (p, e) into (IHTj) where we can again apply Rouche theorem and 
thus find a unique solution P = /3(e). For real e, /3(e), /3 C = /3(e) is the unique solution of 
the system and being self conjugate it is the solution of the monodromy problem and 5 will 
be a rigorous lower bound for the convergence of the perturbative expansion. Obviously if 
we want to optimize the outcome, we have to choose B as to render 5 as large as possible. 
As we shall see in the following choosing a too large B makes the bounds on G increase 
faster than B and thus S has to decrease to satisfy ([S3]) . On the other hand it is obvious 
that B small requires 5 small. 
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As already mentioned the first approach of eq. (1261) requires the integration along a shorter 
path in the z-plane. As the simple bounds on fj and gj we shall give below behave 
exponentially in the length of the integration path, the first approach, even if eq. ( I26j) is 
more complicated than eq. (127)1 . is more apt to give a larger lower bound on the convergence 
radius. 

To compute such lower bound on the convergence radius we use the following simple 
rigorous inequalities 



\fiM - 1| < cosh(v/|e|(|/3|+mi) \co 2 \) 
|flfi(wi) - e(/3oj 2 - CM)\ 



smh(y/\em+ mi ) \u 2 \) 

U 2\ < /, , ,, , == l w 2| 

+ mi) 



< 



co 2 \ [sinh(VNPT^i) H) - ^ 
\g 2 (ui) - 1| < cosh( v /|e|(|/3| + m 1 ) |w 2 |) - 1 



(84) 
(85) 

(86) 
(87) 



where mi is the maximum of the modulus of p(z) along the segment [c<;3,c<;i], and similar 
inequalities for the functions and with argument u 2 . More elaborate inequalities can 
provide a larger lower bound for the convergence radius. We report in Table 1 the lower 
bounds on the convergence radius for a few values of the modulus r obtained with the 
above described method. As expected the square (for which we know that /3 is zero) gives 
the largest lower bound. Due to the exponential behavior of the inequalities ( 1841 - [57)1 the 
bound shrinks to zero in the highly asymmetric configurations. The method applies for 
any r. 



T 


i 


2i 


3i 


U 


5i 




0.1202 


0.05244 


0.02581 


0.01512 


0.00988 


B 


2.1667 


3.80021 


5.64317 


7.435 


9.19403 



Table 1: e c is the rigorous lower bound on the convergence radius. The physical region 
for the coupling is < e < 1/4. B is the parameter appearing in eq.( 



8 ^-perturbation theory 

In this section we shall develop perturbation theory directly from the Liouville equation. 
The analytic nature of the perturbative expansion has to be borrowed form the rigorous 
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treatment of the previous section; on the other hand one can easily obtain in this way the 
value of the action to third order in r\. 

The Green function on the torus of half-periods u% = real, u 2 , r = uj 2 /coi is given by § 



<?(*) = -log 



Vd 



+ 



i(z — z) 2 



16o;i (oo 2 - 0J2) 



satisfying 



AG(z) = 5 2 (z) - - 
a 



*9) 



with a = Aujfrj. The arbitrary additive constant in G has been chosen in (|88p as to have 



J G{z)dz A dz = 



Then we expand 



^>o + 4>i + 4>2 + ■ ■ ■ to have 



1 -AM^ = 2rf V {5 2 (z)-^) 



from which 



(f> x (z) = -8nr]G(z) . 
This is true if the constant term in G is chosen as in (188| 190]) . Next we have 

- tA0 2 (z) = -G(z) 

4 a 



and then 



64ttV 



G(z — z')G{z')dz A dz 1 — . 



(90) 

(91) 
(92) 

(93) 

(94) 
(95) 



Using (I91|93p and 9[(0) = 2rf D , being r/ D the Dedekind modular function, we obtain for 
the X of section 



X = lot 



4w 2 r 7 



477 log 



from which 



— S 2 = 77 — 77 lOg — -J— + 2?7 log 
2,71 ^{Tj 



(96) 



+ 0(ry 3 ). 



(97) 



Eq. (l97|) can be compared with the exact result for the square (15T|) where with e± = — e 2 
1, e 3 = we have 

v^r(i) 



0J\ = —1002 



4r(|) 



(9* 



Here we use for 6*i the convention of [34 not the one of 
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Using 

no(i) = T Si • = ^ (") 

we obtain 

= -77 log 77 + 77 (1 - 2 log S|( - 3 log 2) + r/ 2 log 2 (100) 

which agrees with the expansion to second order of fl5T|) . Similarly one compares eq. (1971) 
with the exact result (155]) for the equianharmonic case finding agreement. 
Using the exact relation 

\dl<P-\{d z <py = e{P + p{z)) (101) 

and the expression of <pi fl93l) one can also retrieve the accessory parameter /3 to first order 
i.e. eq. (17Ti) using 

and the Legendre relation [33] . 

It is easy in this framework to obtain 02(0) i.e. X to second order which integrated 
provides the action to third order in 77 for any r. This would be very tedious to obtain 
in the approach described in section [JJ which however provides the value of the accessory 
parameter to second order. 
We must add to eq. (l96l) </>2(0) given by 

MO) = -^Pr [ G{z')G{-z')dz a dz l - = -^X)'? Wt ~T2 ( 103 ) 

where the prime means m = n = excluded. Such a sum, using standard resummation 
formulas [37], can be rewritten in terms of two simple sums which converge rapidly due 
to the presence of the imaginary part of r 

7T 4 W 1 

S 



E, 1 "K sr-^l 1 
7 W7 =~ ~T? = 77 + 27r / ~T7 — ^(cot7rrn - cot7r™) 
m — rnYim — rnr 45 L — 'nn 6 (T — r) d 

2 , 1 -, 2 (-lt— + -^—). (104) 

f * tj^ T — ti^ Gin 7TTn Gin irm 



If we integrate in 77 according to eq. (l35j) we obtain the third order contribution to the 
one-point function 

Ls? = . (105) 
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9 General analytic properties 



In this section we shall examine the general analytic properties of the accessory parameter 
as a function of the coupling and of the modulus at the non perturbative level. 
We start with the remark that the uniqueness of Picard solution implies also the unique- 
ness of the accessory parameter as 

e^Sge"^ 2 = -Q(u) (106) 

which identifies uniquely 0. Actually can be obtained from 

— / e^dle-^du = € - t <*±l + ^ . (107 ) 

2m J ei 4 (ei - e 2 )(e 3 - ei) 8(ei - e 2 ){e 1 - e 3 ) 

In [SU] it was proven using Green function technique that when e varies in the physical 
interval < e < 1/4, the functions ip, d u <p, d^tp are uniformly bounded functions of u 
in any region of the u plane, obtained by excluding finite disks around the singularities, 
with bounds which depends continuously on e. 

Thus taking the contour of the integral (11071) at a finite distance from e\ we have that is 
a bounded function of e when it varies in the physical region. Such a result combined with 
the uniqueness of the solution implies that is a continuous function of e. In fact if e n is 
a sequence of values converging to eo, due to the boundedness the corresponding sequence 
n must have at least one limit point. However a limit point due to the continuity of the 
basic relations (I18II19P is a solution of the monodromy problem and being such solution 
unique there must be only one limit point. Continuity plays an important role in the 
following as in most of the procedures related to the zeros of analytic functions [3~T] . 
Starting from the relation (1191) we recall that if at a point eo in the physical range 
^12 (C2) / we have also M 2 i{C2) 7^ as explained in section [2J On the other hand if 
Myi(p2) = we have also M 2 i(C 2 ) = 0. We cannot have at the same time Mi 2 {Ci) = 
and Mi 2 (C 2 ) = otherwise the parameter k would be left undetermined against Pi- 
card's uniqueness theorem. Thus given any value €q in the physical region dividing either 
by Mi 2 (C 2 )M 21 (C 2 ) or by Mi 2 (Ci)M 2 i(Ci) we reach in an open interval around e the 
structure 

A(P,e)=B(P,e) (108) 

with A analytic function of and e and B analytic function of and e. 

As done in section [7] it will be useful to employ the polarization technique [36] introducing 

in addition to an other independent complex variable C . 
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We consider now the system 

E 1 = A(P,e)-B(p e ,e) = Q 

E 2 = B((3,e)-A(f3 c ,e)=0. (109) 

We look for solutions of the above system for < e < 1/4. Obviously if ((3,(3 C ) is a 
solution also (/3 C , j3) is a solutions but we shall be particularly interested in self-conjugate 
solutions i.e. those for which (3 C = (3 insofar they are the solution of the monodromy 
problem. Actually from the existence and uniqueness result of the monodromic solution 
we know that for < e < 1/4 there is always one and only one self-conjugate solution, in 
addition, possibly, to other non self-conjugate solutions. In the following we shall denote 
such unique self-conjugate solution as (3(e). 
For e = eo we have 

A(P(e ),e ) = B(Wo),e ) . (110) 
The Weierstrass preparation theorem [351 [3T] can be applied to A(j3,e) — A(f3(e ),e ) if 

A(P,e )-A(P(e ),e ) (111) 

is not identically zero in /3. This can happen only at a finite number of points in the open 
interval otherwise ^ = i.e. A would be a function only of e which from the structure 
of the Mjk(Ci) of section 3 is not true. We exclude such a finite number of points. 
Thus except at most a finite number of points we can apply Weierstrass preparation 
theorem [321 EI] 

A(0, e) - A(/3(e ), e ) = u(/3, e)((/3 - /3(e )) m + c m ^(e)(/3 - /3(e ))— 1 + ... + c (e)) (112) 



S(/3 c ,e)-5(/3(eo),eo) = W (/3 c ,e)((/3 c -/3(eo)r + ^ 1 (e)(/3 c -/3(eo))"- i + ... + ^o(e)) (113) 

with u((3,e),v((3 c ,e) units and Ck,gk analytic functions of e, vanishing at e . 
We consider first the case: m = n = 1. 



At e , /9(e ), /5(e ) we have for the system (11091) the Jacobian 

d(E 1 ,E 2 ) 



J= *,a R V = -H^o),e )\ 2 + \v(/3(e ),e )\' . (114) 

d(P, Pc) 

If J we can apply the implicit function theorem according to which the solution 
and f3 c is unique (and thus self-conjugate for real e) and /3 is an analytic function in an 
open interval around eo and thus we have local analyticity. 
If J = then we look at the equation 



u(re ia + /3(e ), e ) re ia - v{re~™ + /3(e ), e ) re~ ia = 0. (115) 
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For r 7^ divide (11151) by r and call it F(r, a). 



F(r, a) = u(re ta + /3(e ), e )e tQ - v{re~ ia + /3(e ), e )e^ a = . (116) 
Consider a solution a of F(0, a) = 



u(P(e ),€ )e^ =v(P(e ),e )e-* ao (117) 

which is soluble because J = 0. Then in the product of the open intervals ao — 5 < 
a < ao + 5, —5 r < r < 5 r we have that F(r, a) is a C 1 function of a and r, with 
F(0,cto) = and F a (0,ao) ^ 0. Then for small r we have one solution a(r) for a [5%] . 
thus a self-conjugate solution with r ^ (any r in the above interval) in addition to the 
(3 = (3(eo),(3 c = /3(eo)). This however violates Picard's uniqueness result. The conclusion 
is that either J ^ or the Weierstrass polynomials f lll2|lll3[) cannot be both first order. 
In the same way one excludes Weierstrass polynomials (11121) and (I113p with the same 
order, m = n > 1 and \u((3(eo), eo)| = \v(/3(eo), eo)|. 

We go back now to the system fll09j) 
Given eo we have 



J B(/3(eo),eo)-A(/3(eo),e ) = 



A(/3(e ),e )- J B(/3(eo),e ) = (118) 



and in a neighborhood A of eo, of /3(eo), A^c of /3(eo) using Weierstrass preparation 
theorem we can write system (I109P as 



U(/3,/3 c ,e)P 1 (/3 c -/3(e );l3,e) = 



V(0,p c ,e)P 2 (p e - P(eo);P,e) = (119) 
which, as U and V are units, is equivalent to 



P 1 (P c -P(e o );P,e) = 



P 2 (/3 c -/3(eo);/3,e) = 0. (120) 

Necessary and sufficient condition for the two polynomials in (I120p to have a common 
solution in p c is that the resultant [391 ED SO] of the two polynomials Pi and P2 is zero 

R(P 1 ,P 2 ) = f(P,e) = 0. (121) 

In particular we know from the existence result that 

/(/3(e ),e ) = 0. (122) 
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Exploiting again Weierstrass preparation theorem eq. (11211) can be written for e G A x C 
Ao, G A^x c A^ as 

M (/3,e)P(/3-/3(e );e) = (123) 

with 

P{/3 - /3(e ); e) = (/3 - /3(e )) m + ((3- /3(e )) m - 1 a m _i(e) + • • • + a (e). (124) 

In order to apply Weierstrass preparation theorem to /(/3,e) we need that f((3,€o) does 
not vanish identically in /3. The vanishing of /(/3,eo) would mean that the system (11091) 
at e has solution for all /3 near /3(eo). This means, using the Weierstrass-polynomial 
expression for A and B, that m = n and |it(/3(eo), eo)| = |t>(/3(e ), eo)| and this implies the 
existence of infinite self-conjugate solutions with (3 ^ /S(eo) at eo and this goes against the 
uniqueness theorem. 

We start now by computing the resultant R(P, P') i.e. the discriminant of P. If it is not 
identically zero it can vanish in the interval around eo included in the Weierstrass set at 
most at a finite number of points, otherwise it would be identically zero. Thus except at 
those finite number of points we can apply the analytic implicit function theorem [35] to 
have analyticity of /3(e) in a open interval around eo- 

The general case can be treated by computing the reduced Gram determinants D n of the 
power-vectors of the roots [31] 



D r , 



So Si ■ ■■ S n -i 

si s 2 ■■■ s n 

s n-l S n ■ ■ ■ S2n-2 



(125) 



where 



£ + & + ■'■ + & (126) 



being the m roots of P. Being D n a symmetric polynomial of the roots it is a polynomial 
in the coefficients a&(e) and as such an analytic function of e. If R(P, P') = D m vanishes 
identically it means that we have at each e a double or higher order root. Then compute 
D m -i. If it is not identically zero it means that the maximum number of distinct roots is 
m — 1 and the set where they are m — 1 is open and given by subtracting from the initial 
open set the zeros of D m _i. These are isolated points [3T] and thus finite in number. 
Moreover in the region where the maximum number of distinct roots is reached all the 
solutions of P(f3 — /3(eo);e) = (the so called local sheets) are analytic [31j , and in 
particular Picard solution is analytic. 
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Suppose now that D m = D m _i = 0. Then compute D m _ 2 and proceed as above. If D m _ 2 
is not identically zero it means that the maximum number of distinct roots is m — 2; it 
can vanish only at a finite number of point and except at those points all solutions of 
( I12ip are analytic. 

The procedure ends due to the fact that D\ = sq = m. The vanishing of all D n , n > 1 
corresponds to the situation where we have only one m-time degenerate solution i.e. 



from which /3(e) — /3(eo) = — ^a m _i(e) which is analytic. 

Removing the described finite number of points we have that given any eo there is an open 
disk around eo where all the solutions of (j!2ip and in particular the unique self-conjugate 
Picard solution, are analytic except for at most a finite number of points . 
We saw in section [7] that a finite interval around the origin < e < 5 is covered by 
the convergent perturbation theory treatment. For the remainder 5 < e < 1/4 we can 
associate to each e an open set with the above properties and then as 5 < e < 1/4 is 
compact we can extract a finite covering. We conclude that the unique self-conjugate 
Picard solution is analytic on the whole physical region except at most at a finite number 
of points. 

Similarly one treats the dependence of (3 on the modulus r. 
Choose any To belonging to the fundamental region and eo with < eo < 1/4. 
From now on we shall neglect in the notation eo i.e. we shall work at fixed e. 
We start again from the equation 



As done above it is useful to apply the polarization technique to (3 by introducing an other 
independent complex variable (3 C , but this time we apply the polarization technique also 
to the variable r, the modulus, by introducing in addition to r an independent complex 
variable r c . We remark that in the previous treatment of the dependence of (3 on e we 
could have applied the polarization technique also to the variable e but being the physical 
values of e real we would have reached the same results. Here instead the physical values 
of r are in the complex. 
We consider the system 




(127) 



A{P,T)-B(p,f)=0. 



(128) 



A((3,r)-B(f3 c ,r c ) 








(129) 
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We look for solutions of the above system for r in the fundamental region and r c = f. 
Obviously if (3, (3 C is a solution also (3 C , (3 is a solutions but we shall be interested is self- 
conjugate solutions i.e. those for which fi c = j3 insofar they are the solution of the mon- 
odromy problem. Actually from the existence and uniqueness result of the monodromic 
solution we know that for < eo < 1/4 and r c = f there is always one and only one 
self-conjugate solution, in addition, possibly, to other non self-conjugate solutions. In the 
following we shall denote the unique self-conjugate solution as /3(r). 
Chosen tq in the fundamental region we have 

A(f3(r ),r ) = B(M^),f ) . (130) 
Applying the Weierstrass preparation theorem to A((3, r) and B(f3 c , r c ) we have 

A(/3,r)-A(/3(r ),r )= M (/3,r)((/3-/3(r )r + c m _ 1 (r)(/3-/3(r )) m - 1 + ... + c (r)) (131) 

B((3 C , r c ) - B(pfoj, 7b) = v(p e , r c )((/3 c - J{^)) n + ^_ 1 (r c )(^ c - Wo))^ 1 + ... + g (r e )) 

(132) 

with Ck analytic functions of r vanishing at To and gt analytic functions of r c vanishing 
at r an d u and v units. 

As done for the dependence on e, for m = n = 1 if J(tq, f ) ^ we are in the analytic 
situation while J(to, t ) = is excluded by the uniqueness result. In the same way one 
excludes Weierstrass polynomials (I13ip and fl 1 3 2 j) with the same order m = n > 1 and 
K/3(t ),t )| = \v(J}(T ),f Q )\. 
Given r as 

B(^),f )-A(f3(r ),r ) = 

A(Wo), r ) - B(/3(t ), to) = (133) 

in a neighborhood A of r , A c of f , of (3(t ), Ap c of /3(r ) we can write system (11291) 
as 

U((3, (3 C , r, tJP^c - Wo)- (3, r, r c ) = 

V((3, (3 C , r, r c )P 2 (/3 c - ^); /9, r, r c ) = (134) 
which as [/ and K are units is equivalent to 

Pi(Pc-I(^);P,t,t c ) = 

P 2 (/3 c -J(roj;/3,r,r c ) = 0. (135) 
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A common solution of (11351) in /3 C implies (necessary and sufficient condition) the resultant 
of Pi , P2 to be zero 

R(P 1 ,P 2 ) = f(f3,T,r c ) = 0. (136) 
In particular we know from the existence result that 

/(/3(To),T ,fo) = 0. (137) 

/(/3,To,To) cannot be identically zero in for the same reasoning as the one performed 
after eq. (H24j) : thus we can apply Weierstrass preparation theorem and write r, r c ) 
for r G Ai C A, r c G A lc C A c , (3 G A^i C A^ as 

f{/3, r, r c ) = u(P, r, t c )P(/3 - /3(r ); r, r c ) (138) 

with 

P(P - ftTo); r, r c ) = (/3 - /3(r )) m + a^r, r c )(/3 - /^ro))™" 1 + • • • + a (r, r c ) (139) 

and the coefficients a n (r, r c ) analytic in r, r c and vanishing at r , f . 
Thus the equation has become 

P(/3-/3(r );r,r c )=0. (140) 

(3 is analytic in r, r c at all points r, f except those at which P'({3(t);t,t) = 0. These r 
satisfy the discriminant equation 

R(P,P') = D m (r,f) = (141) 

with R(P, P') analytic in r, r c being a polynomial in the a n (r, r c ). 
We distinguish two cases 

1. D m (r,f) is identically zero. Then due to a theorem on polarization [36J D m (r, r c ) is 
identically zero. 

2. Otherwise D m (r, r c ) can vanish only on a t/im set [551 EI], of which the points r such 
that D m (r,f) = are a subset. Thin set have zero measure [35J. Outside such thin set 
the equation is invertible and thus analytic function of r, r c , i.e. /3(r, f) a real-analytic 
function of r. 

In case 1 i.e. 

P m (r,r c )=0 (142) 

we compute D m _i. If it is not identically zero it means that the maximum number of 
distinct roots is m — 1 and the set where they are m — 1 is open and given by subtracting 
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from the initial open set the zeros of D m _x which is a thin set and as such of zero measure. 
In the region where the maximum number of distinct roots is reached all the solutions of 
(11361) (local sheet) are analytic fJT], and in particular Picard solution is analytic. 
Suppose now that 

D m = An-1 = . (143) 

Then we compute D m _ 2 an proceed as above. 

The procedure ends due to the fact that D 1 = m. It corresponds to the situation where 
we have only one m-times degenerate solution i.e. 

P{P - /3(r ); r, r c ) = (/? - /3(r, r c )) m = (144) 

from which we have /3(r, r c ) — (3{tq) = — — a m _i(T,r c ) which is analytic in r, r c and thus 
0(t, f) real-analytic in r. 

We can divide the fundamental region of r in a denumerable set of horizontal strips which 
are compact. We have a zero- measure set of possible non real-analyticity points in each 
strip and the union of such infinite zero measure set has zero measure. 
We conclude that for each e in the physical region the accessory parameter is a real- 
analytic function of r in the whole fundamental region except at most for a zero measure 
set. 

10 Conclusions 

We have considered the problem of accessory parameters on the torus. The specific case 
we dealt with is that of a single source which corresponds to a special cases of the Heun 
equation. We proved that necessary and sufficient condition to obtain monodromy at 
all singularities is the fulfillment of a single implicit equation. Several features of the 
accessory parameter can be extracted from such an equation. A perturbative series was 
developed and rigorous lower bound on the radius of convergence of the perturbative series 
has been given. The second order result for the accessory parameter and third order result 
for the one point function was explicitly computed. 

Modular invariance is useful to find the value of the accessory parameter in some special 
cases and it is satisfied by the perturbative solution. General analytic properties of the 
dependence of the accessory parameter on the source strength and on the modulus have 
been proved. The real-analyticity of the dependence of the accessory parameters on 
the moduli is an essential step in proving Polyakov relation on the sphere which has 
the meaning of determining the response of the on-shell action on the position of the 
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singularities. We shall devote a separate paper to the structure and meaning of the 
Polyakov relation on the torus. 

The described technique can be extended to treat the four-point case or higher number 
of points on the sphere, higher point function on the torus or higher genus surfaces. 
In [18] an integral expression has been given for the accessory parameter for the four 
point function on the sphere. Such a procedure should be extensible to the one point 
function on the torus. However comparison of that result e.g. with the second and third 
order result of sections [3, [8] will not be immediate as our result is an expansion in the 
source strength while the results of [18] are nearer to an expansion in the position of the 
singularities which in the case of the torus correspond to the value of the modulus. 
On the other hand in [201 EE] analytical technique for dealing with the deformation of 
the torus have been developed. These could allow the comparison of the result obtained 
along the lines of [18] . 
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